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abstract 

Four different classes of 

used to eliminate narrow-band int- P 51gnal ca ncelers can be 
^gnal: ( 1 ) cascaded s “conS-oider fr ° m a ^^band 

in-line notch filters, (,t) second-nrd tch filters, , 2) high-order 
and (4) high-order bandpass no“sf c nee, Se can ^lars, 

structure based on second-nrrto . . ^ lers - Of the four, a 
cancel ers is found to perform better^h* f iJ ters used as signal 
The adaptive algorithm for these filfo a °. the other structures. 

>^n and Martin and modify by Pe ^agLL m Pr ° P ° 5ad b * 

/ retragi 1 a , Mitra and Szczupak. 

The Kwan and Martin 

complexity without degrading “I"" be reduced in hardware 
algorithm that out-pe?W a ’ f!”" u * in 0 a ne " « d «Ptiv. 

sa J h nir£d : * F r - 

Shift-Key <BPSK, s:gnals with 



0-0 Introduction 

a" reCe P tlm °* a 

"n*thT^™~ «- 

less ^ergv than the broid-band 5 inn aT'w^ inte ^-^= may have 
concentrated over a narrow bandwidth h ‘ bdCaU5e U 
signal. In the specific wort Lt l^^ masl: out the broad-band 
IS to el i minate only the narrow-band ^ inte \ nttrnt ^ the goal 
Processing will separate broari-hTT 1 nter " fl erence. Later 
signal. However, this later nm P noise "from the desired 
unless all narrow-band ?. 1 1 "° b "rk properly 

broad-band processing is begun. eliminated before the 

interference is to m’ak^u se^f *1 ^scad^? 9 the narrow ' band 
Based on extensive analysis and de of notch filters [1-8] 

ttertin developed a par c .larlv !??* r, r nU1 W ° r,! ’ K " a " and 
modified by Petraalia, hitra !„2 fh “ r ' C5:l whibh waa 
good success. However, the hard! * ,BJ and a PPlied with 
performance with more than two or th requl ^ ed to achieve this 
becomes too complex for easvTmn^ Tx mterfer i ng signals 
high sampling rate is required. One fjh"’ parU CL| 1 a rl y when 
Of this report is to introrim-o , . °* the major contri buti on= 

that maintains the performance of The Kwan ~ Martin algorithm 
requires substantially l e p S h arnV _ h original algorithm but 

band interference when there are mufti pleTnterf erersT 

of narrow-band inte^erencrfrL^ro-d^^d 1 ' 05 ^ 65 l ° eliminat i°n 
describes four different structures ?w signals. Sections 1 
cancellation. Section 2 compass tt^f adaptive interference 
basis for the selection of the second Sactlona and gives the 
the best design. Section 3 gives de ‘ [aT . band ' pa5a canceler as 
algorithm for the second-order band-n^l * htocmati on on the new 
experimental data to demonstrate ?? P canceler and extensive 
algorithms including the algorithm h advanta 9 es °ver other 
dives detailed examples of the new J ^ MarUn - Sedt i°o 4 
Cctncslling narrow-band BPSK di i 9° ri thm s performance in 
BPS f ai gnal s and noi se. Le g 9 ^s'are^ "^^" ° f b - ad - ba nh 

including a discussion of thp TIT Q1Ven ln Sec tion 5 

— Section 4 contain^th^ref erences^ or^the°report ^ eSe 
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1.0 Adaptive Notch Filters 

• „ mui f i ni p narrow-band interference can 

r..~ ..w.- 

:!sr.s^s^ - 

ssr& a- 

Bf&tsrr^as sasr 

Thus IIR pipelining may become an important issue. 



1.1 Second-Order Cascaded Notch Filters 



The second-order notch filter is used 
with the signal as shown in Figure 1-la. 
for the notch filter is given by: 



in cascade and in line 
The transfer function 



1 - 



9— L- 

r-. i 



2(2-k 2 -k 1 ) 



2-k. 



:- x + s ' 2 



H n (z) = 



( 1 - 1 ) 



- (2-k 2 -k 1 2 )z 1 + < l-k 2 ) z 



-2 



For arbitrary values oF H, and k 2 , this is a symmetric notch 
fPter Kith unity gain at DC and the Nyquist Frequency. IF 
kept constant, then the 3db notch width is also kept constant. 
Thus k mav be adapted to remove one narrow-band signal. A 
cascade o^’such inters can be used to remove multiple narrow- 

band signals. 



k„ is 



,2 Second-Order Cascaded Signal Canceler 



1 

The cascaded second-order signal canceler approach shown in 
Figure 1-lb has the advantage that the desired signal does not 
pass throuqh the adaptive Filter. Instead, the band P a *= 
is used to detect the narrow-band signal which is then “btract.. 
From the desired signal, ft constant 3db bandwidth notch can be 
achieved by selecting a band pass Filter wi e 

function: 




<a> CascaiJe ° f Second-Order Notch Filters 




y (n) 



(b) 




Cascade of Second-Order 



Cancel 1 ati on 




Fi 1 ters 



y (n) 





( d) Hiqher-Order q im - = i n 

Signal Cancelling Filt 



Four Categories of Notch Filters 



Figure 1 - 1 . 




k 2 (l - 



(1-2) 



H E P (z ) 



2D (z ) 



where: D(z) is the same 



denominator as H N (z) in equation 



“uSir™ 

sru .» r » ; .« rsuf«5r- 

•frequency parameter k t i -< J * ’ The sensitivity -function 

o-f an adaptive version o-f this -filter. 

H s (z)=H ep (z)*H 55 (z) where: 



H_(z) = 



2k t 2 






- i 



(1-3) 



The par arceter 



k , ma y then be adapted by the formula 



k t (n+D = k t (n) - #e(n)s(n) 



(3) 



Considerable information i. available in_the Hterature^ 
“rihr^di^rhenrU to be best suited tor our 

application. 



3 Higher-Order In-Line Notch Filter 



1.3 

, „ i _ i r ha= the advantaqe that it is 

The notch filter °f Figure 1 ^ be constructed as a 

more easily pipelined. Sue pipeline and has the 

linear-phase FIR filter which FIR filter will require 

advantaqe of linear phase. However, an FIR 
many weights to obtain good performance. 

ftn „ R in-line filter «> * - J ° ^i^ ST" 
rr.^: Pipelining the IIR Filter is .ore difficult than 
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pi pel i ni no thp ftp + 

second-order HR f iiter ^ t “* j“ s than p lpe i ininQ the 

be f^Pted by having a discrete set of ‘ ,™ £ Ulter coild 
could be switched in by a detection • Pr * 5elected notches which 

" ,ax,mums in ‘he spectrum of the input ^ Simply looks tor 
approach could be used with any of th l 9 ?? 1 ' fllth ™9h this 
parti cul arl y attractive with t ■ 'in filters, it is 

generate geometric pole/sero patterns^ h^ 0 ^ * bilit Y to 

' rns m higher-order -filters. 



1.4 



Higher-Order Signal Cancel, 



bene t ^ts^nc^di sadvantacies ^ ^ h “ - the same 

discussed above. However? since th 9 ° rd er. in-1 ine notch 
the signal, the desired = i an - 1 h & cance ^ er is not in-line with 
filter. fl i S o, it J? L'elZr n " d to DaSB through the 

IS -filter. An additional advantaGe^ ^ the adaptive Portion o-f 
adapt the filter off-line and then fr ° m the at >ility to 

15 t0 — «» P ep f ormance oT£ 



1,5 HR vs FIR Realizations 

notch°of L 'bctfon°i . ^and^the 0 second 1 "^ B d enCe ? ^ 5 — d order 

sect on 1.2 would by necessitv T band “ p ass canceler of 

second-order FIR filter cannnf t0 be IIR real i zati ons. a 

band-width in order to adequate?* * SUfficien tly narrow 

1 4 0ret ^ alIy ’ the higher-order reafizati ^ ^ terferin d signals. 
?; 4 Could be either IIR or FIR H Nations of sections 1.3 and 

ideaHy suited for our situation ° Wever ’ the second-order IIR i 5 

resonance with only one parameter the ^ Pr ° vides a ^p 

cequency, to be adapted. Hi a hp r ' ^ nobcb or band-pass 

greac complication without siqnif icant l UR filters would add 

significantly improving performance. 

•; n many appl i cati ons , however if 
IR rather than IIR filters FlR-’fnf^ 7 be advanta Peous to use 
order to realize the narrowband rh ^ ^ be hi 9 h ^ order b 
primary advantaqes of the FTP d ^ haracteri sti cs. The two 
^ ade to be linear phase In* ^ ^ that the V oan be 

Discrete Fourier Transform (DFT } such ^ desi ^ed using the 

FFT) chl P or a Recursive DFT u=inn r ^ FaSt Fourier Transform 

bFT using Residue Number Arithmetic 
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[9]. In the insta ^ ^oul^^designed using one of the 

SnT™ This «** -£V * ^"ome 

linear combiner (ie: tapped-de ay ''"^ype^ o( sectlon 1-4 

more sophisticated adaptive es ■ recursive 

we would use the DFT implemented mth^an^FT chp^ 

DFT to produce a ser ‘ es J? the weights of the output subtracter 
alqorithm would simply adapt the weignus 

that would cancel the various sine waves. 

Figure 1-2 shows an i-P'-^ation of the higher -order FIR^ 

filter of section >- 4 ^ ieh produces N 

FIR filter is set by the initial _ The m reC ursive 

equally spaced zeros arouni d t Qf the de lay ,‘ N ar e used 

pole producing re5ulti ^ in m band-pass filters. 

to cancel filters in turn cancel the narrowband 

These band pa_^ fil adaDtive alqorithm uses 

interference from the input signal. ^Jhe^dap^ive^^g ^ ^ ^ 

a qradient search to find 11 inn all of the narrow-band 

minimize the output power thus cancelling all of the 

interf erence. 
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y (n) 




Figure 



Adaptive Recursi 



ve DFT FIR Notch Filter 
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2.0 Superiority of the Second-Order Bandpass Canceler 



The second-order implementations of section 1.1 p "*’ 

considerable advantage both in ^^Ue^gU-orde^ 
adaptive performance when compa d t t p ed-del ay-line type 

realisations of sections 1. a and 1.4. « W r equine 

in-line notch filter such as that of section i , 0 . weiaht 

Ur 1* adaptive ban^ass 

raUlUsucJ^thk of "^section 1.4 would reguire over 100 
weights to accomplish comparable performance. 



The only higher-order section that appears_to b. c«prt dive 

W ith a second-order U^be required in Figure 1- 

approach is used, a Til ter can be made comparable to 

2, but the adaptive part of the f>‘« performance, the DFT 

that of a second-order system In terms o ^ pa5sband 

approach of Figure 1-* will c ^ eat * . but a5 i on g as N is kept 

ripple in the resulting oupu ^s\ _ ^ j em> Although we believe 

large, this ripple shot. attention our experiments 

that this approach deserves further attention, o P the 

indicate that the second-order band-pass ^ uture 

best choice for our "Irdware ^ developments in such 
developments in real Arithmetic hardware for the 

^“it orRU. !-2 "ight make this approach attractive. 



2.1 Cascaded Second-Order Adaptive Notch Filters 



The filter of section 1.1 was ^the prelU^f 

“ =; rf n r noise" ^uinUTindicated that the cascade 

approach performed VB ^ [ U^^urf 2-1^^^%“^^^ 

used? ^nTnpu? Tonsils of 3 Ue waves at 'U^ne' Uf 
: /T/.0 4=i° and 135°) in Gaussian noise. Each sine wave 

And the noise have the same energy in the Q 

amplitude 1.414 and the noise i 5 Gaussian ^^^^Jiirnotch’ 

fl?^, 2 F?“lf:ho:s P the output of the second stage, ^ 
Figure 2-ld shows the output of the third stage o e 
Instead of eliminating completely one sine wave, the fir 
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w t? c 



iZlfillij-S ines at 36, 45. and 



350 -- 

300 
250 

0 

M 

P 200 
L 
I 

T 150 - 



100 

50 



0 



IjSjegrees in Noise) 



'U-Ar'L 



0 .05 




• 2 -25 .3 3S 7 

frequencv q 



•« .5 



Figure 2-la. 



Input to 2n d-Order 



Cascaded Adaptive Notch 
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M W C= 




Figure 2-lb. 



First Stage of 



2nd-0rder Cascaded Adaptive Notch 
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Figure 2-lc. 



Second Stage of 



nd Order Cascaded Adaptive Notch 
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Figure 2-ld. 



Third Stage o-f 2nd-0rder Cascaded Adaptive Notch 
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° f the *‘-«-Vte x*z 



2.2 The Kwan and Martin Filter 

dete.-ting ^^^^^ 50^31 ^ 153 ’ the pr ° b lem of 

:°uh 2 jeTe:r d r ith a p ~ <* 

z:xTZ' b z f a th b r?*“ ™tp;r i r.^ h \* inu “ ? -* 

Perform both d V this^lr^b^s^uc? 

configuration which is shown to^® h- filter UrB 

extremely well. ° W " to be hl 9hly robust and performs 

the number d of a b. V quad 9 Ject f iorr ^ Martin structure is that 

configuration is given >" thp adaptive filter 

sinusoids to be detected and rein w* N is the nu "’ ber of 

in real time situations with more^Th' I h * S become s impractical 
geometric increase in the reauirert^ 3 "., 4 Slnus ° id s due to the 
p. opose a modification to the Kwan ^ m**" 0 ’ In th * s section, W e 
reduces the required hardware to " U , h Un 5trupturp that ’ 
the performance of the system. ‘ quads with no effect on 




Figure 2-2. 



ne Section of Kwan and Mi 



Notch Fi 1 ter . 



2.3 Kwan and Martin Structure 



The Kwan and Martin structure consists of a “scade.of 

at ace cons^of fb^l- mter liTu “Vh 



1-r, 



l 



1 " 2 



-2 



h b P i ) 



( 2 - 1 ) 



-1 , 2,- 2 

-/ 1 - Zr^osSj z + r i ^ 



where: 



u ; = 



id, = 



pole radius erf the i-th section 

peat frequency of the i-th section 
sampling -frequency 



Kwan and Martin identify two dif f erent ^ethods^f or adapting 

the -filter, ^ ^^^iite^in 1 which the pole radius rj is a 
the “constant bandwidt sdaoted An alternative 

constant and only the frequency u >i _ s ' CUSS ed i n Kwan and 

approach which keeps a constant Q is also disc u 
J H . . pifhpr of these approaches may be used witn me 

structure we will proposing, although too^ave In 

concentrate on the -^^"t bandwidt ^ quan tity in such a 

addition, Kwan and ^rtin Bine the notc h frequency from 

way that it is fairly ea-> ... have no t followed 

the adaptive pa ^™ e ^ e ch " e * 1 the adaptive parameter h-cosS,. 
SublXting ^his in'to equation <2-1) yi.ld. the following: 



1-r j 



1 - 



-2 



•■*E P i ^ 



( 2 - 2 ) 



1 - 2r 4 k i 



^ - 1 



+ r 



2_-2 



r Finnre n— n w e see that the notch f 1 1 ter for each sect i on is 
From Figure ^ «-? we bee <-• , -c-iifor hence: 

the difference between 1 and the ban p -< ’ 



H w j (z ) - 1 _ ki B p i (z ) 



(2-3a) 
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H Ni <Z> 



. 1 - + z -< 



1 - 2r i k i 2 -‘ + r J-* 



(2-3b) 



CALCULATION OF THE GRADIENT 
>he basic structure n-f tk B v 



X(2)0- 




H F PN (c ) 



Y N <2) 



^~|^E F ( N-l ) (Z 



Figure 2-3. K 




and Martin Notch Filter for Tracking N 



ing N sinusoids 
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The overall 



transfer function is given by: 



T(z> 



N 

= n 

i = i 






= n (i p i ) 

i = 1 



(2-4a) 

(2-4b) 






j(z) = L£Az)l 2 = T(z) 2 X(z) 2 



(2-5) 



Hence. , 



the gradient of the objective function J(z) 



is given by: 



6J(z) *T(z> 

= 2E t (z)X(z) 



( 2 - 6 ) 






6 k : 

J 



resps-t to each 1- j : 



61 (z ) 



6k J 



s .in Hi (z) 

i = 1 



i t 



J 



6H nj (z) 



6k j 



(2-7) 



From equation (2-3a) we have: 



From equation (2-2) we have: 



6H n ; <s> 



dk . 

J 



~ ~Hg pi (z) 



-2r :z~ 1 



D ( 



^bp; ( 2 ) H e . (z) 
z) J -J 



(2-9) 



where: 



D(z) 



- 1 - 2r j k jz ~ 1 + r 



: (2 ) = 



J 'j [denominator of H , p ( z > ] 

2r . z - 1 



1 - 2r.k. z -> + r -2 



2 -jil 

D (z ) 



Substituting equation (^-P> ,• + 

. h Lion 9 ) into equation ( 2 - 7 ) 

(z ) 



we obtain: 



ok . 



iJl Hwi (Z) H BPj ( 2>H S (z) 



J i*. 



(2-10) 



‘rte - the comDlete 

would require N-l biouari'- n 1S product for each section 

r- 1 ; »iq U .d. j«t b, t q o a B d ^t^“ r r: ting in * 

e to reduce this by u=ina the n f' r ° duct * Kwan and Martin are 
the input to their cascade J ° utput of the bandpass f ilter ® 

biquads are needed"’^ I °"W (N-J, additiona, J ” 



N 



N 



> t <N- j) = ) N _ 

4 Z I 

J=1 J=i 



/ J = N 2 - 



n 2 +n n 2 -n 



(2-11) 



j = 1 
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Addita this to the N biauads required to realize the cascade of 
notch "-filters and the N bi quads required to realize the H 5i <-> 
lectors Yields a total number of biquads given by. 



$ of bi quads (Kwan/Martin) 



N 2 -N 



N (N+3) 



+ 2N = 



( 2 - 12 ) 



X(z)$ 



2.4 Improved Structure 

Fiqure 2-4 shows the improved adaptive notch filter structure 
proposed in this paper. ^ , x( ,, P for the cascade of the 

^mer^an": ^iUen as' a product 

ir^ii'rtrei^rut'x^r-nrrhro^t^.u, * «- 

bandpass -filters: 




Fiqure 2-4. 



New Notch Filter tor Tracking N sinusoids 
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(2- 13a) 



El<2 ’ ■ i ? l H Ni (*, X(z, 



= X(z) 




(2-13b> 



To get the product of w f \ 

squat: ion <2- 13 b, to .i.pV.« ^Tth^V^ ~ -V 



use 



1 



1 



J x H Ni <Z> X(z) 






T(z) X (z ) - y („ ) 

j 



(2-14) 



X H :'i ^ 

* H BP .( Z; . Thus we have: Stage does not need the 



# of 



bi quads (new structure) = 5 N ~1 



(2-15) 



tFo qU ih ed f ° r the °K«an lnd P Marti’n st nU,,lber of bic I u sd5 

the ZlV?£? t ZZ. H up 
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~n to 






Table 2-1 




cr 



Experimental Results Using the New Structure 



The better^ than^'waf^anc^Mar tin in 

aU ot J the experiments. Figure 2-5 “ ^ 5 %° ne^aves" 

algorithm with 5 notches SPPJ® ws th ; input signal frequency 
imbedded in noise. Figure - u r t r -filterinq and Figure 

response, Figure 2-5b shows the output af ter fi! l ^ t er'-f iltering. 

-5c shows the -Frequency response o Da rameters. Due to 

igure 2-5d shows the adaptation of the five.paramete^^ ^ 

computational complexity, we we ^e un * . being able to implement 
the kwan and Martin algorithm, thus just being ab e / Martin 

the new algorithm represents an improvement over 
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filter °is e uied ZTt^rl LTo'n.V^" th * 5 ~"°^ ^pt.ve 

Fi f^ 2 -^a shows the input spectrum lmbedddd in no.se. 

nutpot, Fi gure 2-6o shows the LtpuLsn'T* 2 ' 6b 5ho “ 5 thd 
shows the adaptation of the 5 notches " 9 tL UI ”’ Snd F ‘ 9ure 2 ~ 6 d 
cancei_the three sines, while 

fi 1 ter "LLsed LeTthtl ^en the 5 -notch adaptive 

2 ; 7 a shows the input spectrum, Figure 2 ° 7 h ' ^ in "° isp - Fi 9 -e 
Figure 2-7c shows the output soprtr 2 7b ^ows the output, 
adaptation of the 5 notches^ What haooe^ 2 ‘ 7d 5hows the 

in eresting and shows just how robust th 5 i""' FlQUre 2 ~ 7d is very 
■first 4 notches lock solidlv nn fK r- th al 9° ri thm is. The 
43°, 60°, and 90°). The -final f G . lrst 4 frequencies (30° 
remaining three frequencies < 20 ° °Lio leF * Jum P in 9 between the 
unable to cancel all three j t L ’ ’ and 150 >• Since it is 

“ch °ne, thus attenuating ’each ^ ‘Sf til “n"' ^ ° f tlme 

scan across the entire range in hope LL . 0dcssl °nal 1 y , it will 
solution. But it quietly setUes h ?V " 9 a better 
attenuating each of the remaining sfnes * Pattern ot 

The -final experiment we triPd with th= 

Figure - e . Here we used 10 notche- and WaVES 15 shown in 

demonstrate the ability ot thTJ ^ 0 Sines in or der to 

-ves. Figure 2-8 shows the La w Pr ° CdSd 10 dldd 

the adaptation of the 10 parameters. 



^•6 Conclusion 



We have demonstrated th^t fK 

structure uses considerable f the new adaptive notch filter 
proposed by Kwan and Martin vet £ ar J? Ware than the structure 
typical applications. I n addi ti onToTh' 5 thG Same ( ° r b0t ter) 
here, chapters 3 and'* 0 ? thi^renc^ d ™° n strated 

Provide other experiments to verffv ttf L° M ’ S ' Th * P15 [1B] 
dbru - bure * This wort has also ho ^ ' P er * orman ce of the new 

C10 -:n. In all Of these c 1 ^ at two conferences 

well or better than the Kwan and hLL i n^sLuct ur e! “ 



i n 



- 21 - 




Figure 2-5a Input of 



5 Sines 



(30 ! 



45°, 36°, 60 



°, and 135°) in Noise 



pmin =0.000001, P max = 100000 0 -' 
= -Cl. 414, 1.414, 1.414, 1.414, 1*V 
ANG = v0.0, 0.0, 0-0i 0 - 0 ’ 



step = 0.100000, s_l imi t = 0.000000, 

Number of sinusoidal inputs is ^ 

Err ...» 
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Input 3 Sines ( 30 °, 45 ° PU 7 6 ° f i C ,o° tch Ada Ptive Filter 
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Input 3 Sines <30° f ^5 %' a S? ° * 5 ~ Not< = h Filter 
j ep = 0.100000, s lifl|i ’and 13^ ) m Noise 

S"S, »■« ■ ■««»«>.»» 

“ rif‘ "iif '.-•■• . ,:i» 




Fiaure 2-6a 



Input o-f 3 Sines (^0°, 36 ? 



and 135°) in Noise 



step - 
Number 

FREQ = 
DC = 0.000? 



- °-~- p 

of r vrr 1 \~- 

NOISE = 1.000, k_l imi t ~ 1 ^ 00 _’ 



pmax 
1.4143 
0.0003 
1, n 



= 1000GI 



sec = 



i s m = /a 000 0.000, 0.000, 0- 000 ^ 0 - 000} 

r e a dii - ' F «- S!5«. «■*». 

Second order -forgetting. Factor - ,0.^, «. 



0.9003 
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Figure 2-6c Frequency Response o-f Filtered Output of 5 Notch Filter 



Input 3 Sines (30°, 36°, and 135°) in Noise 



step = 0.100000, s_limit = 0.000000, pmin =0.000001, pmax - 100000. 

Number of sinusoidal inputs is 3: AMP = {1.414, «’ rarapn. 

FREQ = C0.083, 0.125, 0.375} ANG = {0.000, B.OMU 0 - 000 ' _ , 

nr = Pi 000 NOISE = 1.000, k_limit = 1.000, slow - 1, n_sec 

Initial values of kl = {0.000, 0.000, 0-000, 0 - 000 ’ 0 - 000> 
i ^ • • _ /g O00 0 900, 0.900, 0.900, 0.900J 

Seconder der% or get ti ng. Vector = <0.900, 0.900, 0.900, 0.900, 0.900) 
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pmin =0.000001, pmax = 100002 



step = 0.100000, s_l i mi t = 0.000000, 

Number o-f sinusoidal inputs is 7. 

AMP = 11.414, 1.414, 1.414, 1.414, 

FREQ = (0.083, 0.125, 0.167, 0.250, 

ANG = {0.000, 0.000, 0.000, 0.000, 

DC = 0.000, NOISE = 1.000, k_l i mi t = 

Initial values o-f kl = C0.000, 0.000, 

Notch pole radii = {0.900, 0.900, 0.900, 0.900, 0.900' 

Second orcer forgetting. Factor = C0.900, 0.900, 0.900, 0.900, 
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Output 5-notch Martin Filter (7 sine input in noise) 

— - — i r— r” 1 — i 1 

300 



250 




Fiqure 2-7c Frequency Response o-f Filtered Output o-f 5 Notch Filtei 
incut 7 Sines (30% 45”, 60° , 90°, 120°, 135°, and 150°. in Noise 



pmin =0.000001, pmax = 100000- 



step = 0.100033, s_limit = 0.000000, 

Number of sinusoidal inputs is 7. 

AMP = <1.414, 1.414, 1.414, 1.414, 1.414, 1.414, 1.4 4, 

FRED = f 0 . 08 T - 0.125, 0.167, 0.250, 0.3oo, 0.C.75, 0<4 ^ J 
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0.900, 0.90. 
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Fiaure 2-B Input o-f 10 Sines in Noise 
(30 ., 36% 45% 60% 72% 108% 120% 135% 144% 150 ) 

step = 0.100000. ..U^t - 0.000000, pmin =0.000001, pmax-10W»« 

"IT. -"1/ ^3 >•-; i:£; fcSfc 2! 

TS : £». VZ\ S:U; S: «; y-. T\ 0 ~: ? 

DC = 0.000, NOISE = 0.000, U" 11 ' l - 000 ’ 

e!m, 0.000, 0.000, 0 - 000 - 0 - 000 ’ 0 ’ 000 ’ 0 ‘ 000 ’ ^ 

NOt f 0 .C O?900, = 0.900, 0.900, 0.900, 0.900, 0.900, 0.900, 0-900, ».«• 

0^0:900, 0-900, 0 .900, 0*900, 0-900, 0-90) 



3.0 Detailed 



Description o f the New Algorithm 



new algorithm and demonstrate it= performs^ SCU5S detai 1 the 
The derivation of the structure and L 5f ° n simulate * data, 
discussed in sections 2.3 and *? 4 u* Qradlent calculation were 
details -for implementing the adant*f H 5ha11 give the 

choose the various adaptive otrZtt ’ lnformat *°" on how to 

effects of different choices on the perfo^ dem ° nstrati of the 

tne P er f ormance of the filters. 



3.1 The Adaptation Algorithm 



struc^n^c^JirLoXpero^s: ^ ^ aI ™m. 

types of second-order filters: 



The 



Band-pass filters ( 



see equation 2-2): 



H tc; <z) = 



1 - z " 2 

2 rr^ r . l; . 



(3-1) 



ar, d sensitivity filt 



ers (see equation 2-9): 



H c • (z) = 



2r . z ~ 1 



1 - 2r. I:.--* + 



- 1 
. z 



D (z ) 



(3-2) 



where: 



D(z) 



= 1 - 2r . k . z ~ 1 
j j 



+ r . 
j 



t denominator of w 

"BP i *2 / J 



The partial derivative is given by (see equation 2 10) 



d T(z> 




H EPj (z)H 5j (z) 

j 



(3-3 



) 



The only parameter that is modified in the adaptation process 
is the pol e cosine k- Tor each bandpass -filter section j. >nis 
parameter also effects the denominator of the sensitivity filter 
H . (z). The only other parameter , r-, is fixed. Selection o , 
will be covered in section 3.3. The standard LMS algorithm would 
give us the following update procedure for adapting k^: 



dT(z) 



k : ' = k : - lie - (z) 

j j j 



(3-4) 



dk j 



where 



: a is a fixed step size (typically 4-0.1) 



As derived in chapter 2, the correct choice for the error 
function e(z) is: 



= (z) = T(z) (for all j) 

_ j 

However, an alternative choice which can be derived by forming 
the gradient of the section outputs Yj(z) is: 



e j ( z) = Y j (z) 



(3-6) 



Both of these choices for e^z) were tested and found to work 
wel 1 . 



A word of 
to imply that 
sense , we are. 



explanation is in order here. Equation (•> 4) seem 
we are mixing z-domain and time domain terms. In 
The derivative are taken on the frequency domain 
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actuany U tiIe e ser 1 t es U <Z l /dk i in Fi iur e 2-4 are 

s^HiBSr r-rti- 



3.2 Forgetting Factor 

5 d ofu\"o°n 3. . 

x ol f to test vari ous forqetUnq f actor/" 0 "" ^ 5ee L52) ‘ lr > 
following update equation: ' 1 s ’ we employed the 



f-J HO : -k j Ck J-step*gradi ent* ( 1 . 0+ p n ; ; } 



0 J/pmax)/ (pmi n+p Cf; , 0 ] ) 



where: 



step 


* (fixed for a particular structure) 
dT(z) 


(3- 7a) 


gr adi ent 


= e . (z) 






V 

ok . 

j 


(3-7b) 


FLk,Q] 


output of forgetting filter 


(3-7c) 


pm ax 


- limits the effect of large pCI , )0] 


(3-7d) 


pmi n 


= limits the effect of small p n : , B: 


(3-7e) 



NOTE: 



the 



^ -Found that pmax and 
Lpdate algorithm could 



Pmin were not required 
be replaced by: 



so 



in reality 
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(3-8) 



k j [ k ] : =k j C k ] -st ep*gr adi ent/ptk ,03 



However, all of the data in this chapter was run with equation 7 
setting pmaK=1000000 and pmi n=0. 000001 . 



Three different forgetting filters were tested. The first, 
called zero-order forgetting, uses the sum of the squares of the 
i nstantaneous outputs iKz./ik, as pCk.M. The next -Ued 
first order forgetting, is a simple sing e-pole first order low 
lUter aDDlied to the section output Ci>T U ) /hU. ■ 1 . I 

third type, called second-order forgetting, is a second-order 
low-pass filter applied to CAT (<- ) /dk ■ 3 . 



Zero Order Forgetting 



d r k , 0 3 




dT(z) 




(3-9a) 



First Order Forgetting 



p [ k. , 0 3 = (l-r f ) 



dT(z) 



dk, 



+ r f p £ k , 0 3 



(3-9b ) 



Second Order Forgetting 



p l k , 0 1 



output of one of the following second-order filters 
with input CdT(z)/dkj3 : 
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simple 2nd-0rder 



Low— Pass 



H t (z) = __ 



h i <1 ~ 2s 



- 1 



, -2. 



1 - 2 r l - “ 1 

* V 



+ 2 



-2 



<3-9c) 



Complex 2nd-0rder Low-Pass 



H 2 (z) 



where: 



h 2 (1 ~ 2cos (55 j ) 2 -1 + 2 - 2 j 
1 - 2r.k -2" 1 + -.'2 

f j - 

e j = ArcCos (k • ) 

J 



(3-9d) 



The filter of equation r?-o-n 

3-3 Selection of Notch Bandwidth by and Mai "tin C 5 ]. 

Each bandpass filter h« - x • , 

parameter r If we let » Z \f 1Xed band width set by the 
notch and Ae the" £ W 3db f ^-ncy ot the 

bandiadth is d = x y . wfU al ? ? nCy °* th * notch - then the 
frequency 2 =j ( 90 o f .. =0) ~ ^ flr jt calculate d at the 

other values of k/baied^n « y for 

from, the i in owfng'equaU on: ' " ^ tW ° 3db * ^quenci es *. and y. 



- = 1 H n (c ) P = 

O 1 


" 1 +r 2 ' 


2 


U+cos2et) 2 + si n 2 2wt ] 




o 




. <l + r 2 cos2 u t > 2 + r’ sin * 2 ut_ 



<3-10a) 



1 = 


'l + 2 r 2 + r * ' 




4 + 4 cos 2 ut 




4 

. 




. + 2 r 2 cos 2 ot + r * 
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(3-ii) 



cross-multiplying yields: 



l+ 2 r 2 cos 2 ut+r 



f _ 



= i+: 



lr"+r 



* + cos2ot (l+2r 2 +r 4 > 



Using the fact that cos2ut 
-for cosut: 



= 2cos 2 ut-l, we solve equation <3-U> 



-Ok 



cos2ut = 2COS 



! ut - 1 - 



(3-12a) 



1+r 



1-r' 



(3- 12b) 



cosot 



- + 



2 (1+r 4 ) 



Since x e is the lower 
■frequency, we have: 



3db -frequency and y c is the upper 



v = Arc Cos 
c 



2 v * 



-(1-r") 



-j 2 (1+r 4 ) 



(3- 13a) 



v_ = ArcCos 
' 0 



r i -^ 2 



-j 2 ( l+r 4 : 



(3-13b) 



A1 so 



, 51 nce sin 2 e + cos 2 0 = 1 , we have 



1-r' 



cos(Xg) - " 



2 ( 1+r 4 ) 



(3-14a) 
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(3- 14b) 



cos(y ) = + 



1-r' 



-[3 



(1+r 4 ) 



5in(;.; 0 ) = + 



1+r" 



3 



2 (1+r 4 ) 



si n (y 0 ) = + 



l+r‘ 



33 



(1+r 4 ) 



(3- 14c) 



<3-14d) 



^ "' a/ h “ CaICUl3te ^ bandwi dth d fro . d . „ . y . 

C > 0 ■ 



cos (d) = cos(x e -y .) 



= cos(x e )cos(y B ) + sin (>: c )=in (y 0 ) 



(3-15) 



Substitution from equations ( 3 - : 



14) int ° e 3 u ati on (3-15) yields: 

cos (d) = ~ (1 ~ 2 r 2 +r') + (l+2rV 4 ) 



2r 2 



2 ( 1+r 4 ) 



1+r' 



(3-16) 



the * . notch , t 
"•%> < “ s ; ny v,:zv n h i 
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z" 2 = 0 



( 3— 17a) 




2rk 



1 - 







= 0 



(3-17b) 



Substituting -for z=e J f =cos (fr ) + jsi n </7 ) 
setting the real and the imaginary part 



in equation (3- 17b) 
to zero yields: 



and 



4rk ; 

1 - 005 (/t) + 005(2/7) 

1+r 2 



0 



(3-lBa) 



— sin (/r) + sin ( 2 / 7 ) - 0 



Usinq tri qonometri c identities: 

4rk • 

l - - COS(/7) + 2C05 2 </7) -1=0 

1+r 2 



(3-18b) 



(3- 19a) 



4rk . 



i sin(/7) + 2sin (/7)cos(/7) - 0 



1+r' 



(3- 19b) 



both equations yield the same result: 



cos (fr ) 



2rk 



1+r 2 



(3-20) 
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pole X«cV “lightly different froo the 

va ue of r. From the *«»■•.* of d^o the ‘ 

»ell suited to the type of notrhff ! vaUles of r around .9 ere 
bend interference. The teble also n!t eliminate "arrow- 

the uctual notch frequency / th<= i 1 the measure d values for 
upper 3db frequency y for^ee’h of “Ho^h * , and £ 

olosely to the calculeted veluef^L * SS \ These c ™P^e 
^ ’ anc * as follows: ua lon ^ 16), equation ( 3 - 



Table 3-1 Effect of r 



on Notch Width 
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hence: 



y 2 - xy -dy = y 2 - f 2 -dy 
y = d/2 + (d/2) 2 + / 2 



0 



x = d/2 



(3-22a) 

-j (d/2) 2 + / 2 (3-22b) 



3.4 Minimum and Maximum Frequency of Notch 

Equation (3-20) implies that the notch frequency is slightly 
different from the pole frequency B. In particular. 



2rk j 2r 

cos ifr) = = 7 CD5(S) 

1+r 2 1 +r* 



(3-20) 



This difference becomes smaller as r approaches one. 
the pole frequency 6 is at its minimum of 0 , then k 



However , 
=1 and: 



if 



\ = f 7 



- Arc Cos 
0 = 0 



2r 



1+r' 



If 6 



is at its maximum of it (180°), then kj 1 and. 



(3-21 a) 



tt - A = /V 



= ArcCos 
0=TT 



-2r 



1+r' 



(3-21b) 



From equation (3-20) we see that the difference bet th 

frequency fr and the pole frequency 6 is at most A. “ 

approaches s/2 (90°), the difference between the notch and pel 



frequencies goes to zero, 
various values of r. 



Table 



II gives £ and 180°-£ f or 



Table 3-1 1 Effect of r on A 




dre lm P°rtant because if there is a f qUencies Q^en in Table 3-1 1 
thc,n the minium notch frequency or ^ n ^° ld at a frequency lower 
notch will be forced to thi= minimum. 15 a DC ^lue, the 

there is a sinusoid above th- ' frequency. Similarly, if 

the " 0tch wiiI be%or^d' J to ° r at the 

Flgure 3 -> =bo W s the result of l DC ! " laximum frequency, 

process. The notch is forced to 6 0 -^^ . ° n the Captation 
the raaius of the notch filter is rhnt 10 thiS CS£e beca ^e 
Figure 3-2 shows a signal ^ thp m ? a5 r=0 * 9 ‘ Similarly, 
notch to the maximum frequency of forcin g the 



3.5 Algorism Performance 



Critical to the perfo-manrp ^ 4 -i~ 
the various parameters: ' hS al 9°nthm is the choice of 
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1. Number of notch -filters, 

2. Notch width r 



4. 

5. 

6 . 



Notcn Wl QUI r j , 

Zero, -first, or second-order forgetting (or none), 
Step size fi, 

Limit on the parameter k, and 
Limit on the parameter p. 



The first three choices are somewhat independent of each other. 
But the choice of the step size a is highly dependent on the 
®ir, t three choices. We found that limiting k such that - .^1 

was' very important in all cases. However, we found that limits 

on p were not very useful. 

The choice of the number of notch filters has a profound 
effect on the algorithm performance, but is assume ° * ° t 

our control. We assume that the number of notches is fi>>ed at 
rme number by the application. Thus in evaluating the 
algorithm, we look at three cases. 



Case 1: Performance with only one notch, 
Case 2: Performance with five notches, 



P -s 



8.50 o 



Performance with ten notches. 



T,bl. ’.-Ill summarizes experiments with the new algorithm. 
For thesJ experiments, the following 10 frequencies were used. 



8 

9 

10 



f 

f 

f 

f 

f 

f 

f 

f 

f 

f 



i n 
i n 



l n 
i n 
i n 
i n 
i n 
i n 
i n 
i n 



= 0. 08333 f 
= 0. i f 5 
= 8.125 f 5 
= 0. 1667 f 5 
= 0.2 f s 
= 0.25 f s 
= 0.3 f 5 
= 0.3333 f s 
= 0.375 f 5 
= 0.4 f. 



(30°) 

(36°) 

(45°) 

(60°) 

(72°) 

(90°) 

(108°) 

( 120 °) 

(135°) 

(144°) 



-n cases where only one frequency was selected, the T® su ' ts 
are based on the results with each of the 5 (1, 3,^6 

8. 13) taken separately. In cases where three freque 

used, the following five combinations were use t , 
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s'KL-Ht " 



O a 



Cone! usi ons 



fh •-■ IV shows the results of fh 

he-e results, we conclude that a tornff f ■ 0Xperi ments - From 
multi-notch filters and that ! for 9 ettlr >g factor i s needed 
does not seem to inprove the perto°" d ° rder for getting factor 
factor which seems to work u Zl T ^ * fir ^-order 

order^f or get ti ng-f actor is the ** The first- 

si tu P t • C ° nvergenc:e - Here we refe^to^ d ° eS ° 0t e>:hlbit 

la=^ not n h WHen tHere are more interferersTh COnver 9 er >ce as the 
— notch m the cascade "jumne- ho f ^ han notches and the 
frequencies to be eliminated Th • J; ween the remaininq 
actually be desireable in some in=ta JU, " P " conver 9enc e ; miqht 
-t to avo,d the „ nst-ocdXXui “* ^ ““ - -U 

detLn Se i al ° ne ’ th e e p^Imeters U p the a ndT° rder f ° r9ettin 9- ^or 

r. 

and 

convergence for the zero-order *' f? llmits to obtain 

--tr ba 

z J ---ssit :;Xt 
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Table 3-IV Summary o-f Tests on New Algorithm 




1 . 



In all cases without noise, parameters converge to value ^ 

o-f notches is greater or equal to number of sines. Iterations listedi 
table represent convergence to three decimal places o corre 

In all cases with noise, parameters oscillate around exact value.^Tb 
number o-f iterations listed in the table is e nu 
required to establish this oscillatory pattern. 
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p «in =10 6 and 



^th er types 

3 .ax=10* 6 >. 



- 49 - 



4.0 Algorithm Performance and Simulation 

The adaptive digital filter algorithm as described in the 
ear 1 ie r chapters is simulated and tested using synthetic data. 
Simulation was carried out on VAX 11/785 computer using Fortran 
77 The synthetic data generated for testing this algorithm is 

o-f -four categories: 



1. Sinusoidal signals with white gaussian noise 

2. Narrow Band Noise with white gaussian noise 

3. Bi -Phase Shift Keying (BPSK) sequence 

4. Frequency Shift Keying (FSK) sequence 



4.1 Sinusoidal signals 



For qeneratinq sinusoidal signals P la ^ d at different 
normal izid frequencies 8 , a second order AR process given as. 



*„ = 2cos(9 i )K , R _ : 



k-2 



(4-1) 



was used with poles on the unit circle. This approach was chosen 
inorder to reduce the computational burden. Using equation 41, 
we find that initial conditions are very important and they are 
chossn such that x ,=0 and x_„=-sin <« i > giving a unit amplitude 
sinusoidal signal. The 8 , value is between 0 to 180 ano n is the 

number o 



Siqricil. me u , 

frequencies desired. The required signal y k needed to , 



l iumu-i w i — 

input into the adaptive algorithm is given as: 



j = n 



Yk 




i = l 



(4-2) 



where T k 



is a white gaussian 



noise 



N (0 , a 2 ) . 
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4.2 



Narrow Band Noise 



The narrow band 
equation: 



si gnal 



is generated using 



the di f f erence 






= 2rcos (8 



) 2 ' 



k-1 



r 2 - j 

r A k-2 + 



U 



( 4 - 3 ) 



where 6 ■ dec i Hoc 

deSlred S1 9nal y k is obtained w a: dlffsr ^ > "stents!^ The 



i = n 




i = i 



( 4 - 4 ) 



4.3 



Bi Phase Shift Keying Sequence 



The generation of BF'Sh' <=i n r-i u 

are generated as shown in Parts “ hlch 

Tol i ows: gure 4 1- The three parts are as 



by passing a uniformly di ^buted^"" <RBS ’ achie ^<i 

limi ter. (fin important note is thatch" * hard 

the two consecutive bits of RBS is !//) * Val betwee " 

Generation of another sequence ■ 
art in practice a spreading - -/ blnar Y numbers which 
The specific sequence usecMnV a r SRreadlr,g sequence. 
system is normally not avail ahtJl communi ca ti ons 
i ntended receiver' 8ny ° ne but tha 

have generated the spreadinn ' 1Cular 51m ulation we 
uniformly distributed noi s/th by passin 9 a 

hit interval is j/y _ ) rough a hard limiter. The 

Generation of the nh^ Q 

given binary signal which is^he'e* 9 ; * hB mappin 9 of ‘be 

<2) ab °^ into either 0 or w ^ aa„T '! ^ ° f U) a " d 

at appropriate sample times. 
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DSSS BPSK Generator 



Uniform 

Random 

Number 











05 



RBS 



X(i) 




Signal 



X(i> 



(l 


• kf ) 






' 0.5 




Fig. 4.1 



Fiqure 4-1. 



Generation of the BPSK Signal 




<see T Ffg°re P 4-lf hard is stored in 

the array y w PUt of the sec °^ hard limH ‘ * 

the array ^deJ ofy e n ZT-!T 'h * a subt e llTT, " 

ss^^HSS* M* » * ■•”** 

sr^ " «?&aar - tnxrs sr* 

the equations: 1Dn ' The ^ed signal no„ is g - v ‘„ by 



wP k - AcoB(2ir/ e k + « k ) 

= C>: (i ) ® y ( j) j n 



P = n 




P=l 



(4-5) 



(4-6) 



(4-7) 



where i and j are indices n-r th 

•■S an important point note thatTh^ ** defi ? ed earlier. It 
defined by the parameter p = <)** ^ completely 



C » 



/c. / 



b y • 



additional feature of sprang tff ® BPSK 51 ^ n£l but it has an 
^ . 1 p frequency and carrier -frt=n & spectrum by controlling the 
diagram of the scheme f givef.T'f ^ baud ^te. The bfcf 
Vk is given by: - 9^en ln Figure 4-1. The desired sfgnal 



= 



cc . 



+ ^ 



Where: 



1 ^ ^be narrow band BPSK 



(4-8) 



Parta of the frequency speotr, ^h 5 * 9 "* 1 Pl*ced at dtfferent 

^ generated^ I p*"^* ^ brMd »“ d ^ 



4.4 Frequency 



Shift Keying Sequence 



. . n fhic cpnuence needs a random binary intelligence 

distributed ntuse throug a har such that i=k/ b where 

5 ; 0, ;f t e baud-rate" of ^the" in-format ion and k is discrete sample 
'ntUblr Now the desired s.gnal is generated via: 



5 k = 2cOS(S K )S k _ 1 S k-2 

0 k = 8 + £x (i > 

Vk = 5 k + y x 



(4-9) 

(4-10) 

(4-11) 



where 6 is the carrier -frequency 
-frequency modulation. Initial c 
they are chosen such that s_ t -0 
amplitude sinusoidal signal. 



and 6 is the depth of the 
onditions are very important and 
and s_ 2 =-sin(0) giving an unit 



4.5 Simulation 

Ft t;F STi'Si;™ trtt-ttt.rr";™. .. 

Fortran 77. Furthe, deta u5Ed are the Kwan and 

^t-nm^n/r^ned Kwan -d Hart.n tilter 
algorithm) as described in chapters m and The fey P 

of these filters are: 



(a) 

(b) 

(c) 

(d) 

(e) 



Sharpness of the notch filter defined by pole 

Btep^size'in ^he'i^rementaUon procedure (a) 

Time constant of the fading filter 

Model order (n) tnfnrfprers (m 

Order of incoming signal or number of interferes 
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^•5. 1 Fading Filter 



In order to t r r I* f k D +.* 

to compute <l/N)T\ e 2 ' over '^ _Varyin 9 Parameters it is essential 



= \ \/ 

A * k - 1 



+ < 1 



- X)i 



(4-12) 



^/■p^tt/sjVacto^or m ° 5tly ref erred to as the 

to the first-order forgett’n^uat^-^ 5 *’ 

iizt'Xfj!" - - < 0 , by 

optimizing u and X. The value nf achieved by manually 

alUB ° f “ was chosen as 0.01 and Xs0.9. 

using^Narrow Band^Noi se' signal y a & X . t thB al9 ° rithm «» tested 
signal was used only to tune t hi" ? Performance. This 

the notch filter) . It wL ob erve'd thaf f °* 

close unity variance o-f the Da ™f+ th *i lf the value °-f r is 

ased on simulation the value of r was h BStlmate increases. 

ue ° + r was chosen as 0.9. 



4.5. 



Check VAX version against IBM-PC version 



to the^ajrr:: y * was appi - b 

s: fi s: vi 

shows the results of a sine-wave test f 2 ~ 4 ’ Fi 9“re ^ 

chapter 3. The purpose of this test f *’ la5 described in 
version of the proqram oave the same WS5 I? verify that the VAX 
version used in the previous sect lts as the IBM-PC 

consisting of three sine wa^es in the 'ore' ^ 4 ' 3a WaS the in P ut 

th t f' 9nals - Figure 4 -3b is the filto ? enCe °* noi5e i hut no 
that the filter did remove the sfne output, demonstrating 

parameter adaptation. By mahino a d t f’ Fi9Ure 4 - 3c sb ows the 

th^i! ° n b ° th the Wx ^nd the S IBM-Pc ailed comparison of this 
that the VAX program was worlino the ’ ^ WerS abIe *o verify 

sorting the same as the IBM-PC program. 






New Structure 




Figure 

(NOTE: 



4 - 0 . Block Diagram o-f New Algorithm. 
Same’ as Figure 2-4, except tor notation) 



- 56 - 








Figure 4-3. 



Three Pure Sinusoids in 



Gausian Noise 
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4.5.3 Elimination of narrow-band BPSK interference 

Our next experiment was .^ r ^han sine 

algo'ithm performed on narrow input consisting of 

wave interferes. Figure * ^ ^ ^ < requencie5 as the 

three narrow ban • - p . 4 _ 4b 5hows the -filtered output 

sine waves in I Figure eliminate narrow-band 

verifying that the new 4 _ 4c sho „ s the parameter 

adaptation o this case and should be compared to Figure 4-3c. 



4.5.4 Transient behavior (tracking a moving interfered 

The transient behavior of the adaptive filter was tested by 

. S 5 i anal using the hardware of Figure 4-1. The 

generating an FSh sign . b t the spe ctrum is 

FSK signal is similar o e ^ ^ ' spectrum switched back 

constantly changing. In our example, l ^ show5 the FF T of 

s d . ^ ssr ^ 

5how= tne oucpu- ■ 4 _-, c shows the parameter 

two t r thr^i^:^s input 

^ •H«in.Jog .t extremely well C1M. 



4. 



El imi 



ination of narrow-band BPSK from broad-band BPSK 



Our next experiment was aimed at eliminating narrow-band 
BPSK. interference from a broad-band BPSK s1 ^ - 19 =i ls 

shows the input consisting of three narrow-band BFSK.i goals 
interfering with a broad-band ^signal ^^/^""ergy 

content? US Figure°4-6b shows the output ^—^^h^he^ 

signal C1B1. Furthermore, we processed th l* *'!nd obtained the 
Cyclic Spectral Analysis Software Package E v . and 

output shown in Figure 4-7. This outpu was as x p ect-»d 
confirmed that the new adaptive notch filters P 

well C121. 
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- Q -U t P 1 1 f fiPPCt ri|m 




4 00 



»f r n 

»:«%•: ?““« r 3 ‘- 

Xi.-j:;. 1 ;: ° f - 

IdI\ a ? Ce ° f the noise’ '"°'i 5 o 
Adaptive filter step 0 ’5l 

Adaptive filter r ' ' ' J'J 1 




Figure 4-4. 



Spectrum of 



3 Narrow-Band BF'SK 



signals in Noise 
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200 300 400 




input FSK signal 

Samples per carrier cycle of the 

mark and space .... 6 ‘ , „ 

Variance of the noise •••••• 

Adaptive filter step 
Adaptive filter r 0 • * 



Figure 4-5. Tracking of an FSK Signal in Noise 
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I 

0 . 5 Ilf' 'WMW^'y^V 



1000 2000 2 0 00 ,000 



Input . 

Number of ' interf erers " ' ' I ? S ~ BPSK 
interferes r C * rrl « <** «e 

interf erers ^ of '^ °-‘' 5 '^ 

c»°rr d erXu 9nal ' Sa " P1 -- 200 ' 200 

Adaptive filter step q ' 02 

Adaptive filter r "°'° 2 



Figure 



^ 6. Filtering 

BPSK signal 



Narrow-Band BPSK 
in white Gausian 



•from Broad-Band 

Noi se 
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3.82M -i 




a. 



Broad-Band BPSK without Interference or Noise 



11. SHI -| 



e ee: 




b - 



Broad-Band BPSK Corrupted with Interference and Noise 
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4.5.6 



Performance with too many notches 



there were too many notches for^h" WaS . what Would happen if 
signals. Would the remaining notrh numi * er oi interfering 
road-band BPSK signal, and if so eIiminat e the 

Spectral Analysis Software Package F e * fect the C ^lic 

thas important issue. In Figure ts tlT** 4 " 8 4 ~ 9 ad ^ess 

and three notches and in Figure 4-9 t h = tW ° inter terers 

and three notches. As before Fin h ? re 15 onl V one interferer 
spectrum, Figure 4-8b shows the o'*^ 4 8a shows the input 
s OW5 the parameter adaptation. As^an^e^™ 1 "’ and Fi ° ure 'l-Sc 
the remaining notch wanders around thf -t Seen . ln Fi Sure 4-8c, 

° e of the broad-band BPSK signal att rac ? u encies of the main 
However, since the notch is V J V at tempting to cancel it. 

oi the b road-band BPSK signal th* ^ Compared to the main lobe 
any significant amount of^icina^ 15 Unable to eliminate 

the Cyclic Spectral Analyst I WeT ^ 
results as Figure 4-7 [ 101 . a similL^ WE QDt the same 

4 9. The two remainino niches*!?^ reSuft KSS Found in F-'gure 
main 1 obe of the broadband BPsVsia'TVt^ — d the 

eliminate any significant amount ' 9 1 ’ but are unable to 

Cyclic Spectral Analysis s“ wV Pa V,T Th u « *■" the 

agein got the results of Figure 4-7 ri0?f “ appIied - we °"ce 



4.6 Conclusions 
'he BPSK signal 



end the resul ts’showed that 9 the"new the new al 9 0rit hm 

eble to perform well under varTnnf ? PtlVe n °tch filters were 
real world. The new algorithm S situations expected in the 
mission and was able to'prepar'e^th 2154 *" 1 - Y Fared wsU in »ts 
y the Cyclic Spectral Analysis Sof twar^Packag^ ^ prPcessin 9 



given as input to the new algorithm 
- the npw . ... y uim 





1 . 5 






0 . 5 




V-w/^v/ 



Model mismatch m-2 and n ":L_ BpSK 

Input 

Number of interferers ••• 

Samples per carrier cycle o 

interferes I2.0,3.u 

Sampl.es per Chip bit o 200 

interferers 14 

Broad band signal samples per 

carrier cycle 4 • 

Broad band signal samples per 

chip bit • 

Adaptive filter step 0.02 

Adaptive filter r 0.9 



1000 2000 3000 4000 



Fiaure 4-8. 



Two Interferers and Three Notches 





— S,C£j^tXiUIL 



100 



200 300 



400 



Model .ismatch ..j and n . 3 
?“■“>« of l 0 terfe rers ';"“ S ' BPSK 
interfered •' 

interfered Chip ij ^ °i 

cdded^id oarsa " pi «- 200 ' 200 

0Mp d bit° d Si5 ° al sa "P‘es 4 per 

Adaptive filter step . . ! 0 02 
Adaptive filter r [ q 9 



One Interferer and Three Notches. 



- 65 - 



Figure 4-9. 




5.0 Conclusions 



i-^-rs?=.r= asa. »— 

conclusions: 



1 . 



Of the many possible algorithms available for solving 
this problem, only two seem to meet the requirements 
and aUo feasible in terms of hardware complexity 
and hardware speed. These two approaches are. 

adaptive second-order notch filters based on a 
modified Kwan-Martin approach or 



a. 



b . 



an 



approach based on the DFT. 



Ci i 



T hi = report has studied carefully the adaptive notch 
filter approach based on a modification of Kwan and 
Martin and found it to be extremely well suited to the 

problem. 

Single interfered are eliminated within at most a 
-few hundred iterations, 

with multiple interferers, the first notch will 
eliminate one of the interferers with n a few 
hundred iterations and other notches hick in 
quickly such that as many as ten interfer 
be eliminated in as few as 500 iterations, 

c When there are fewer notches than interferers, 
notches can be designed using proper forgetting 
filters to either eliminate the maximum number of 
interferers or to attempt to reduce all 
interferers while eliminating some of them, 

When there are more notches than interferers^ the 
extra notches do not adversely effect the broad 
band signal, 
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e. 



f . 



abffto handU^requency 1 Z^thltTre T 

fnterfer^f^/fh eIi ” inate two or more 

togetllr aV SB lnterfererS « doee 

IBM-PC^r^^fi^com 5 f SS1Iy implemented on and 
hardware implemented oil ^ ^ 5 “ eH 5uited *°r 



5. 1 



Brief Look at Hardware Feasibility 



Although hardware desinn i = k _, 
we shall look briefly at oossihi the 5COpe oi thls report 

order to safety ourselye= that tte'f""" ""Rations in' ’ 
desired application. Ideally we w ls fea5lb le for the 

notch filter that would op^e It T , 1 '* to de 51 9" - 3 to 10 
MHc. Hi th currently available terh 109 * rec l ue ncies up to 100 
available commerciaily todav w™ ih 1 ° ? 9y ’ about the best 
follows, we shall attLpI to iZlt ^ t0 20 In what 

discuss how throuqh the use of pi nan • ’ S 10 20 MHz figure and 

sampling rate to as much as eight tiI«"?i- WB h COUld increasa this 
ue: 80-160MHZ). 9 tlmes the basic sampling rate 



5-1.1 Commercially Available Hardware 

MHz sampIIIg'IItrilluI^^II^rdware feasibility of a 10-20 
mters ’ Ke 5ha11 >«* - tive^r-, 1 ;;;— - b f e ‘- 0 c notch 

!• MI P S R 30 l 0 . 40ns cycle, 2 cycle add - 

and 19 cycle divide C14], ’ ° Cycle mult iply, 
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3- T«» Instruments ^sg-e ’ 

cycle multiply, and 11 cycle aivi 

root ) C 141 , 

i rurlP add. 2 cycle multiply, 

4 . AT&T DSF'32c, 20ns cycle, *- cycle a , 

and 3 cycle divide [15], and 

5. Honeywell HD* 46110, f-^J^^TansTuniofy' 
product of 8 elements for an effective 

and accumulate time [16]. 

The first four chips are of sim *j^ k a ^^ 

could be implemented using 6 ^"^and would require us to 
chip is really an array proces ’ ta ae of its speed. Although 

vectorize the algorithm to advantage * ^ ^ i5 beyond 

the b ;cor^ h this h r:port, so‘we shall concentrate our attention 
*or the time being on the other tour chips. 

Figure 5 -, shows the bloci: diagra^tor^ implement ing^a ^ 

bt ock^ diagram^ makes use" of °a^coet f i c i ent memor y wh i ch can be ^ set 

s£H sss rr±r, 

can then be repeated for all sec c;_ 7i . Jn the process of 

regular structure as indicate in 9 add - tional advantage of 
designing this, it was realised t •- t-wan-Martin algorithm is 

the modified algorithm over the origin,! ^ 

that the original algorithm couid no cascade chain in 

identical blocks because of the need for the 
calculating the filter derivative [ 10 ]. 



1.2 Time Budget 

The details of the time budget have been worked out inthe 



syswrs. Tsjtni s^.nsa'ss’ws 

DSP-32c which yields a sampling rate o - , 

sufficient to meet our 10-20 MHz design requirement. 



This is not 

However , 
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Figure 5-1. 



Block Diagram of 2nd Order IIR Filter 
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Basic Building Block 




6 Output to noxt otogt 

Fig 5.2 




Figure 5-2. 



Block Diagram of 



Basic Building Block 
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3-Notch Canceller 




Fig 5.3 



ure 5-3. 



Block Diagram for 3-Notch Cancel 



er 
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this was done without pipelining the IIR ^‘ion. Usjng the 

, . . fnrhniniips o-f Soderstrand and Loomis LI/ -1J, we wouiu 

ounl^his saving rate by at least two end probably tour 
resulting in the desired 10-20 MHz design. Furthermore, if we 
were to vectorize the algorithm, the Honeywell chip cou 
ootentially yield 100 MHz sampling with both vectorization and 
oipeHmng! Based on this hardware analysis, we feel confident 
that sampling rates of 25 MHz or more can be achieved in practice 
with commercially available DSP chips or most certainly with a 
custom VLSI design. 



5.2 Future Work 



There are two avenues o-f -future wor 
before making a hardware commitment -for 



k that should be pursued 
the narrow-band -filter 



project: 



1 . 



DFT based approaches to the problem should 
especially in view o-f the availability o-f 
FFT processors and such exotic approaches 
DFT's usinq QRNS arithmetic C9], and 



be studied, 
high-speed 
at recursive 



, detailed investigation o-f hardware implementations of 
he modified Kwan-Martin algorithm should be made 
ncluding simulation and layout of a custom VLSI 
ipproach Because of tbs regularity of the structure 
n Figures 5-1 though 5-3, the custom VLSI may well 
yield some very excellent filters. 



1 
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